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Abstract 

We tested the hypothesis that density-dependent competition influences the 
evolution of offspring size. We studied two populations of the least killifish 
{Heterandria formosa) that differ dramatically in population density; these pop- 
ulations are genetically differentiated for offspring size, and females from both 
populations produce larger offspring when they experience higher social densi- 
ties. To look at the influences of population of origin and relative body size on 
competitive ability, we held females from the high-density population at two 
different densities to create large and small offspring with the same genetic 
background. We measured the competitive ability of those offspring in meso- 
cosms that contained either pure or mixed population treatments at either high 
or low density. High density increased competition, which was most evident in 
greatly reduced individual growth rates. Larger offspring from the high-density 
population significantly delayed the onset of maturity of fish from the 
low-density population. From our results, we infer that competitive conditions 
in nature have contributed to the evolution of genetically based interpopulation 
differences in offspring size as well as plasticity in offspring size in response to 
conspecific density. 
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Introduction 

Propagule size is a key life-history trait that directly influ- 
ences both offspring and maternal fitness. Under many 
conditions, larger offspring have higher individual fitness 
(Bagenal 1969; Tessier and Consolati 1989; Fox and Cze- 
sak 2000; Janzen and Warner 2009; Marshall and Keough 
2009; Jacobs and Sherrard 2010; Rius et al. 2010), but 
there are some circumstances in which the benefit of lar- 
ger size is either very small or undetectable (e.g., Kaplan 
1992; Fox et al. 1997; Bashey 2006). As offspring size is 
often inversely related to offspring number, variation in 
offspring size wOl also create variation in parental fitness 
through the tradeoff between offspring size and number. 
High variation in offspring size among conspecific popu- 
lations and related species (Reznick and Endler 1982; 
Amat et al. 2001; Marshall and Keough 2007; Reed et al. 
2009; Riesch et al. 2009; Collin and Salazar 2010) raises 



the questions of how much of this variation is directly 
adaptive, under which circumstances might it be adaptive, 
and for whom, parent or offspring, does the bulk of any 
benefit accrue. 

Evolutionary theories that consider offspring size 
explore conditions under which females should partition 
resources between offspring size and number to maximize 
parental fitness (Smith and FretweU 1974; Brockelman 
1975; Lloyd 1987; Winkler and Wallin 1987; Kindsvater 
et al. 2010). Theory predicts that the optimal parental 
investment per offspring should change when the rela- 
tionship between offspring size and offspring fitness dif- 
fers across ecological settings (Parker and Begon 1986; 
Sibly and Calow 1986; McGinley et al. 1987). Population 
density is one ecological factor that theory predicts will 
play an important role in selecting for variation in 
offspring size and number (Brockelman 1975; Parker 
and Begon 1986). In high-density populations, when 
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competition is intense, selection favors the production of 
larger, competitively superior offspring (reviewed in 
Stearns 1976), especially when the effects of density affect 
juveniles more than adults (Sibley and Calow 1983). This 
prediction assumes that larger body size at offspring inde- 
pendence confers higher fitness at higher densities. 

Empirical evidence supports the idea that population 
density can affect the relationship between offspring size 
and fitness. In plants, large seeds and seedlings have higher 
fitness in high-density environments (e.g.. Harper 1977; 
Stanton 1984; Winn and Miller 1995). Fewer studies have 
examined this relationship in animals, but their results 
uphold this general pattern (Berven and Chadra 1988; Per- 
rin 1989; Hutchings 1991; Einum and Fleming 1999; Mar- 
shall et al. 2006; Bashey 2008). Maternal adjustments of 
offspring size in response to conspecific density (or indica- 
tors of density) also suggest that density is an important 
agent of selection on offspring size. For example, when 
food is limited, females often produce fewer but larger off- 
spring (Gliwicz and Guisande 1992; Reznick and Yang 
1993; Reznick et al. 1996; Bashey 2006). Studies that 
manipulate adult densities (Kawecki 1995; AUen et al. 
2008; Leips et al. 2009) and chemicals associated with vari- 
ation in density (Burns 1995) find that larger offspring are 
produced when females experience high-density condi- 
tions. An adaptive interpretation of this pattern is that the 
larger offspring are better competitors in the higher densi- 
ties that they will face. Despite the attention given to off- 
spring size variation and the adaptive significance of this 
variation, few studies test both facets of this issue (but see 
Reznick 1982; Reznick et al. 1990, 2001; Bashey 2006). 

Here, we examine the relationships among offspring 
size, density, and components of offspring fitness using 
the least killifish {H. formosa) from two natural popula- 
tions in north Florida, Trout Pond (TP) and the Wacissa 
River (WR). These populations differ dramatically in con- 
specific density (Leips and Travis 1999; Richardson et al. 
2006), with population density in the Wacissa River often 
exceeding that of Trout Pond by as much as sevenfold 
during the breeding season. Fish from these populations 
also produce offspring that differ substantially in size; 
wild caught Wacissa River females have fewer but larger 



(>40%) offspring per clutch (Leips and Travis 1999; Sch- 
rader and Travis 2005) and these differences have a 
genetic basis (Leips et al. 2000). However, size-specific 
allocation of mass devoted to reproduction, after adjust- 
ing for differences in female body size, does not differ 
between populations. Thus, the differences in the num- 
bers and sizes of offspring produced reflect evolved differ- 
ences in the way that reproductive mass is packaged, 
suggesting that these traits have responded to population- 
specific selection pressures. 

Our study addresses three general questions related to 
the adaptive significance of offspring size variation in dif- 
ferent competitive environments. First, do offspring from 
high-density populations perform better (exhibit higher 
growth rates, mature earlier, and/or mature at larger body 
size) in high-density conditions than offspring from 
low-density populations? Second, do offspring from 
populations that differ in density differ in their density- 
dependent competitive ability when placed in competition 
with each other? Third, how important is the relative 
difference in body size in determining the outcome of 
competition? When females from both of these popula- 
tions experience high adult density, they produce larger 
offspring than females held at a lower adult density (Leips 
et al. 2009). We took advantage of this plastic response 
of offspring size to social density to manipulate the size of 
offspring from the Wacissa River population. This allowed 
us to examine the role of relative size versus population of 
origin in determining the outcome of density-dependent 
competition. It also allowed us to evaluate the extent to 
which maternal adjustments of offspring size is an adaptive 
response to variation in competitive environments. 

Methods 

Study organism and populations 

Heterandria formosa is a live-bearing topminnow in the 
Poeciliid family that tends to aggregate in aquatic vegeta- 
tion in shallow water. Their small size and external mor- 
phological features indicating sexual maturity make 
H. formosa an excellent species for life-history studies. In 



Figure 1. (A) Adult female Heterandria 
formosa. (B) Adult male H. formosa. 
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females (Fig. lA), sexual maturity is visually indicated by 
the appearance of a small black dot on the anal fin 
(Fraser and Renton 1940); they mature between 8 and 
10 mm (standard length) in approximately 30-40 days. 
In males (Fig. IB), sexual maturity is indicated by the 
development of the gonopodium, a modified anal fin 
used for sperm transfer (Constanz 1989); they take longer 
to mature (40-50 days) and mature at a larger body size 
(10-14 mm) than females. 

Parental fish for this experiment were obtained from 
the Wacissa River (WR) in Jefferson County and Trout 
Pond (TP) in Leon County. These populations are 
approximately 50 km apart and genetically isolated (Baer 
1998; Soucy and Travis 2003), although molecular data 
suggest that this isolation is fairly recent (Baer 1998). 

Animal husbandry 

We measured offspring performance traits using second 
generation (F2), lab-reared offspring. Parental stocks of 
WR and TP fish were initiated with 50 wild caught 
females and 20 wild caught males from each population. 
Breeding was carried out in 19L aquaria in a temperature 
(31°C) and light (14:10 day:night) controlled laboratory. 
Parental densities were maintained at 5-6 females and 
2-3 males per tank during breeding of the parental gener- 
ation, and tanks were checked weekly for offspring. First 
generation (Fi) offspring were removed and assigned to 
an Fi density treatment (explained below). Pedigrees were 
kept for all offspring so that inbreeding in the FjS could 
be minimized. Only individuals that were born in differ- 
ent tanks were allowed to mate. 

Females in the WR FI breeding tanks were kept at one 
of two densities in 19L aquaria (either 6 females and 2 
males per tank, or 2 females and 2 males per tank). These 
different maternal densities are within the range of natural 
densities experienced by WR fish (Leips and Travis 1999) 
and influence offspring size of both WR and TP females 
(Leips et al. 2009). All TP FI breeding tanks had 6 females 
and 2 males, thus producing large TP offspring. Fish were 
fed ground Tetramin® (Tetra, Melle, Germany) flake food 
twice daily ad libitum, so any effects on offspring pheno- 
type from different female densities should have resulted 
from density differences alone, not from reduced per capita 
resources. To summarize, three types of F2 offspring were 
produced: TP offspring from tanks with high female den- 
sity, WR offspring from tanks with high (WR[H]) density, 
and WR offspring from tanks with low (WR[L]) density. 

Experimental design 

Eight treatments were used in a two factor, fractional fac- 
torial design (Dunson and Travis 1991; Warner et al. 



1991) (Table SI). Treatments evaluated offspring perfor- 
mance (growth rate, age, and body size at maturity) at 
either low or high density (8 and 32 individuals) in treat- 
ments containing either offspring from a single popula- 
tion (TP, or WR(H)) or offspring from both populations 
in competition (WR[H] with TP, or WR[L] with TP). 
Thus, there were four single population treatments and 
four mixed population treatments. Single population 
treatments evaluated performance of the TP and WR(H) 
offspring at low and high densities. Two of the mixed 
population treatments measured the performance of the 
larger WR(H) offspring when in competition against the 
smaller TP offspring at both densities. The other two 
mixed population treatments compared the performance 
of WR(L) offspring when in competition with TP off- 
spring at both densities. In this case, the offspring were 
more similar in body size. In each mixed population 
treatment, TP and WR offspring were mixed 1:1 such that 
there were either 8 or 32 individuals per unit, matching 
the total density of the single population treatments. 

Execution of the competition experiment 

Prior to the initiation of a set of experimental replicates, 
all F2 offspring were removed from Fi breeding tanks. We 
then collected offspring for assignment to the different 
treatment groups for the next 8 days. Thus, the offspring 
within each treatment were similar, but not identical, in 
age. A drawback of this procedure is that it minimized our 
ability to detect age-specific trait differences among treat- 
ments because trait measurements were taken weekly. 
However, our results should not be biased because all 
treatments used offspring within this range of ages. A sub- 
set of offspring from each treatment was photographed to 
determine initial body size. In the mixed population treat- 
ments, we tagged fish from one of the populations with 
calcein to identify the origin of the fish in subsequent sam- 
pling events. WR fish were marked in three of the repli- 
cates, TP fish in the remaining two. Calcein binds to 
calcium in tissues and can be seen in fin rays and scales 
(Leips et al. 2001). H. formosa were tagged by maintaining 
individuals for 24 h in a solution of 250 mg/L solution of 
calcein, buffered to a pH of 5.3 using tris (hydroxymethyl) 
amino-methane. This dye is retained for at least 5 weeks 
in the laboratory and had no adverse effects on survival or 
growth compared to an equal number of unmarked con- 
trol fish (Leips et al. 2001). 

One set of the 8 treatments was started in aquaria in 
the laboratory at one time, and was then moved to the 
experimental mesh bags placed in large (850L) cattle 
watering tanks in an open field at the Florida State 
University Greenhouse Facility. Each tank contained a 
mixture of well water and rainwater as well as aquatic 
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vegetation and plankton from Trout Pond. Tank communi- 
ties were established in five cattle tanks 8 months prior 
to the initiation of the experiment. Approximately, 
2000 mL of suspended plankton from Trout Pond were 
added to each tank monthly during the course of the 
experiment. 

Each cattle tank contained a set of 8 mesh bags in which 
the fish from one replicate of each of the 8 treatments were 
raised. Bags were made from nylon tulle material (1.5 mm 
mesh size) and framed with glass rods to maintain a rect- 
angular shape. Bags were hung from cords attached to a 
square PVC frame, and suspended in the tank so that each 
contained approximately 16L of water (Fig SI). Each bag 
contained a small amount of submerged vegetation, M. 
laxum and H. verticillata, to provide cover and substrate 
for epiphytic growth. H. formosa is omnivorous, and fish 
in the bags presumably fed on plankton as well as the veg- 
etation and algae in the bags. 

Each treatment was replicated 5 times (one replicate 
per cattle tank). Each time a set of replicates was started, 
all 8 treatments were initiated in separate 19L aquaria (in 
16L of water) in the laboratory at the same time. For the 
first week, fish in the low-density treatments (cells 1,3,4, 
and 7 in Table Al, Table SI) were fed twice daily with 
10 mg of flake food (an ad libitum amount of food). 
Fish were fed 15 mg twice daily for the second week, and 
then 20 mg if individuals were retained in the lab for a 
third week (see rationale below). High-density treatments 
(cells 2, 5, 6, and 8 in Table Al, Table SI) received the 
same ration per aquarium (i.e., one quarter of the food 
ration per individual). This confounds differences in per 
capita food amounts and density, but was meant to 
reflect the effect of increasing density in natural popula- 
tions. Thus, both resource availability and differences in 
social environment may have influenced the measured 
traits. 

Fish were moved to a cattle tank when they had grown 
too large to escape through the mesh bags (2-3 weeks). 
Each replicate was then censused every week. This 
sampling interval reduces our power to detect subtle 
differences among the treatments but was necessary for 
logistical reasons. Males, females, and immature individu- 
als were sorted by visual inspection and taken to the lab. 
Calcein tagged and untagged fish were sorted by examin- 
ing anesthetized individuals under an epifluorescence 
microscope (Leips et al. 2001). Fish were anesthetized in 
a low dose of MS-222. In the third sampling period, all 
calcein-tagged fish were retagged with calcein to ensure 
that the tag was retained until the next sampling period. 
Fish in the different categories (tagged and untagged; 
mature male, mature female, and immature individuals) 
from all treatments were photographed separately in 
water-filled bins containing a metric scale and then 
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returned to the cattle tanks. The size of fish in each group 
was measured from photographs using an image analysis 
system. 

The entire experiment was carried out over the course 
of a single summer The monitoring period of each repli- 
cate ran from the date the replicate was set up until the 
first female in each bag had reached sexual maturity (or 
the first TP and WR females had matured in the mixed 
treatment bags). 

Statistical analysis 

The initial sizes of the offspring from all three treatments 
(TP: n = 29; WR[H]: n = 19, and WR[L]: n = 10) were 
compared using analysis of variance, with population as 
the categorical variable (Zar 2009). Tukey's HSD proce- 
dure was used for post hoc comparisons of means among 
the groups. We measured the initial size of fish only in 
the first replicate to minimize handling of newborn 
offspring. As the effect of maternal density on offspring 
size matched that of our previous experiments (Leips 
et al. 2009), we assumed the maternal density treatment 
had similar effects on offspring size in the remaining 
replicates. 

We assayed the effects of the different treatments on 
offspring performance using three traits: the average 
growth rates of individuals in each group for the first 
2 weeks, and the age and size of the first maturing indi- 
vidual of each group from a replicate. Growth rates were 
calculated only for the first 2 weeks because growth rates 
of males and females begin to diverge after that time 
(Parrino, Leips and Travis, unpubl. ms.). Focusing only 
on growth rate for the first 2 weeks minimizes the effects 
of chance differences in sex ratio among replicates that 
could have obscured differences among treatment groups. 
Average trait values of maturing individuals from each 
replicate were not used for analyses of age and size at 
maturity because there was large variation in the propor- 
tion of individuals that matured in each replicate. This 
confounded age-structure with density effects among 
replicates. In addition, in some cases (particularly in 
high-density treatments), maturation occurred over a long 
period of time, skewing the mean and median values 
differently in the different replicates. 

We carried out four analyses of variance for each trait 
using contrast codes to explore the role of maternal den- 
sity, offspring density, and population of origin on 
offspring traits. AH analyses (designated as A-D) and 
associated contrasts are detailed in Table SI. Within any 
single analysis (A, B, C, or D), there are three contrasts i, 
ii, and iii. The first two contrasts (i and ii) test for main 
effects and the third contrast (iii) tests for the interaction 
between two main effects. 
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Analysis A uses cells 1, 2, 4, and 6 (Table Al, Table 
SI): the treatments with offspring from only one stock 
population (either TP or WR(H)), at two densities, low 
(N = 8 individuals) or high {4N = 32 individuals). For 
any dependent variable (growth rate, age, or size at matu- 
rity), contrast i tests the effect of quadrupling the density. 
Contrast ii tests whether TP offspring differ from WR 
(H). Contrast iii tests whether the effect of quadrupling 
density differs between stocks, thus testing if the popula- 
tions are differentially sensitive to the depressant effect of 
density. 

Analysis B uses cells 1, 2, 3, and 5 and only uses data 
from TP offspring. Contrast i tests whether TP traits dif- 
fer when they are competing only with other TP fish ver- 
sus when half the fish are WR(H), given a constant 
density. Contrast ii tests whether there is an effect of total 
fish density on TP traits. Contrast iii tests whether the 
effect of density on TP traits depends on the identity of 
the competitor with whom the TP fish are competing. 

Analysis C asks the same questions as analysis B, but 
instead focuses on WR(H) offspring and provides com- 
plementary inference. 

Analysis D uses cells 3, 5, 7, and 8. For this analysis, 
we focus only on data from the WR fish when in compe- 
tition with TP fish. Contrast i tests whether WR(H) and 
WR(L) differ in performance in the presence of TP fish. 
Contrast ii tests the effect of increased density on WR fish 
performance. Contrast iii tests whether the density effect 
differs between WR(H) and WR(L). 

We did not compare the competitive ability of TP off- 
spring when pitted against WR(L) and WR(H) offspring. 
This is because TP fish differed from WR(L) fish in two 
ways, by population of origin and by the density of the 
maternal environment. Thus, population of origin and 
maternal environment would be confounded in this com- 
parison. 

A drawback of this analysis design is the number of 
statistical tests on the same dependent variable; however, 
this approach has some virtues in that light; within any 
single analysis (A, B, C, or D), each contrast is indepen- 
dent using Type III SS in the two-way ANOVA. 

AU statistical analyses (described above) were carried 
out in SAS ® software, V 9.2, SAS Institute Inc. (Cary, 
NC) (ANOVA). All dependent variables were log trans- 
formed to satisfy assumptions of ANOVA. 

We checked for differences in survival among treat- 
ments using logistic regression (Proc GENMOD (SAS), 
binomial distribution, logit link, and scaled deviance 
options). We found no significant differences in survival 
among treatments (x^ = 10.11, P = 0.18) so we did not 
consider this further. The data on age and size at matu- 
rity presented below are only from females because very 
few males matured in any treatment. 



Results 

Initial offspring size 

The WR(H) offspring had the highest initial body sizes. 
TP offspring were the smallest, whereas the WR(L) off- 
spring size was intermediate between the two (Fig 2). 
These differences in offspring size were significant, 
(F2,55 = 7.72, P < 0.01), with WR(H) significantly larger 
than TP, and WR(L) statistically indistinguishable from 
the others. 

Analysis A: the effect of density and 
population origin in single stock treatments 

Quadrupling the density of offspring significantly reduced 
offspring growth rate by about 50% (Fijg = 7.41, 
P = 0.009, Fig S2) and this effect was similar for both 
populations (interaction effects: Fi is = 0.01, P = 0.93). 
There was not a significant interaction between popula- 
tion of origin and density (Fijg = 0.01, P = 0.93). 

The age at maturity for females did not appear to differ 
between densities (Fug = 0.33, P = 0.57, Fig S3), with 
population of origin (Fi_i5 = 0.54, P = 0.47), or their 
interaction (Fijg = 0.30, P = 0.58). Ages at maturity were 
well within those estimated from otolith samples from 
our field sites and independent laboratory studies 
(unpubl. data). 

Increased density had no significant effect on the size 
at maturity of females (Fj^ig = 0.84, P = 0.36, Fig S4). 
Across densities, fish from W(H) stock were not signifi- 
cantly different for this trait (Fi_i6 = 1-12, P = 0.30), nor 
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TP WR(L) WR(H) 

Treatment 

Figure 2. Initial standard lengths of juveniles measured prior to 
assignment to a particular treatment. TP = Trout Pond offspring, WR 
(L) = offspring of Wacissa River females held at low density, WR 
(H) = offspring from Wacissa River females held at high density. 
Shown are average values ± 1 SE. 
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was the interaction between population and density 
(Fi,i6 = 0.15, P = 0.70). 

Analysis B: the effect of density and the 
identity of the competitor on Trout Pond 
offspring 

The first contrast in this analysis was designed to test 
whether the identity of the competitor had an influence 
on TP traits. As expected, increased density significantly 
reduced growth rate (Fijg = 13.68, P = 0.0007, Fig S5). 
We did not detect an effect of the identity of the compet- 
itor on growth rate (TP vs. WR(H) fish, Fijg = 0.24, 
P = 0.63) nor did we find that the effect of density 
depended on the competitor type (TP vs. WR(H)), 
Fi,i8 = 0.66, P = 0.42). 

Quadrupling the density did not affect the age at matu- 
rity of TP fish (Fi,i6 = 0.45, P = 0.51). However, com- 
pared with the values in pure TP treatments, the presence 
of W(H) fish caused the TP fish to delay maturity by 
about 24-33% (Fj.jg = 6.10, P = 0.02, Fig. 3). This sug- 
gests a competitive asymmetry between TP and WR(H) 
fish, with the WR(H) fish being the stronger competitors. 
There was no evidence for an interaction between density 
and whether TP fish were competing with members of 
their own or the WR(H) population. 

For size at maturity, none of the contrasts were signifi- 
cant (pure vs. mixed stocks: Fijg = 0.03, P = 0.86; den- 
sity: Fi_i6 = 0.27, P = 0.60; interaction between density 
and identity of the competitor: Fijg = 0.00, P = 0.90). 
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Figure 3. Trout Pond (TP) females delayed maturity winen tlney were 
in competition with relatively large Wacissa River (WR(H)) offspring 
compared to when they were raised with only other TP offspring 
(comparisons are averaged across densities). Shown are average 
values + 1 SE. 



Analysis C: the effect of density and the 
identity of the competitor on Wacissa River 
offspring 

This analysis was designed to test whether the identity of 
the competitor had an influence on any of the traits of 
WR(H) fish. As in all previous analyses, increased density 
substantially reduced WR(H) offspring growth rate 
(Fi_i9 = 7.25, P = 0.01, Table 1, Fig S6). The identity of 
the competitor had no significant influence on WR(H) 
growth rate (Fijg = 0.05, P = 0.81), and there was no 
evidence for a density x competitor (WR(H) vs. TP fish) 
interaction (Fj^ig = 0.01, P = 0.90). 

The patterns in age at maturity of WR(H) fish when 
alone or with TP fish were slightly different than those 
seen in TP fish. First, increased density actually decreased 
age at maturity as opposed to the lack of an effect of den- 
sity seen in the TP fish (Fig. 4). Second, they matured 
only slightly (and nonsignificantly) later in the presence 
of TP fish, compared with the pronounced delay in TP 
fish in the presence of WR(H) fish (compare values in 
Table 1, Figs 3 and 4). None of the contrasts was statisti- 
cally significant (pure vs. mixed stocks: Fi 17 = 0.10, P = 
0.76; density: Fj 17 = 1.48, P = 0.23; interaction between 
density and identity of the competitor: Fijy = 0.0, 
P = 0.98). These results are consistent with WR(H) fish 
being stronger competitors than TP fish. 



Table 1. Means and standard errors of traits of Wacissa River (H) off- 
spring (from Wacissa River females held at high maternal density) and 
Wacissa River (L) offspring when in competition with Trout Pond off- 
spring (means for treatment groups compared in Analysis D). 

Treatment 





WR(H) 


WR(H) 


WR(L) 


WR(L) 




versus 


versus 


versus 


versus TP 




TP Low 


TP High 


TP Low 


High 


Trait 


Density 


Density 


Density 


Density 


Growth 


0.10 (0.01) 


0.06 (0.01) 


0.11 (0.01) 


0.07 (0.01) 


Rate 










(mm/day) 










Age at 


35(4) 


31 (1) 


39 (5) 


34(4) 


Maturity 










(days) 










Size at 


10.3 (0.8) 


8.8 (0.5) 


9.6 (0.5) 


*9.9 (0.7) 


Maturity 










(mm) 











*The average value shown here contains data from one replicate that 
may be an outlier. The average size at maturity for this treatment 
group without this replicate is 9.2 mm. 

TP = Offspring of Trout Pond Females held at high maternal density. 
Data shown are average values (-1- 1 SE) 



) 2013 The Authors. Published by Blackwell Publishing Ltd. 



953 



Adaptive Differentiation in Offspring Size 



J. Leips et al. 



40 



38 



IS 
T3 



36 



S 34 



32 



ns 

0) 30 
< 

28 



26 



WR(H) alone 
WR(H)withTP 



Low 



High 



Density 



Figure 4. Mean age at maturity of WR(H) offspring (± 1 SE) in single 
stock versus mixed population treatments (means for treatment 
groups compared in Analysis C). TP = Offspring of Trout Pond 
Females held at inigin maternal density, WR(H) = Offspring of Wacissa 
River females held at high maternal density. 



The patterns in size at maturity for WR(H) fish 
reflected those seen in the TP fish, in which none of the 
effects were statistically significant (effect of competitor: 
Fi_i7 = 0.01, P = 0.91, density: Fi^y = 3.13, P = 0.08, 
interaction: Fi_i7 = 0.28, P = 0.60, Table 1, Fig S7). 



Analysis D: the effect of relative size of 
Wacissa River fish on competitive ability 

Analysis D focused on traits of Wacissa River offspring 
reared from the two maternal environments (WR[H] and 
WR[L]) and placed in competition with Trout Pond fish. 
This analysis was designed to reveal how different initial 
sizes and prenatal experiences of Wacissa River offspring 
influence growth, age and size at maturity in competition 
against a common "tester" stock (TP) at low and high 
density. 

The average values for all traits are shown in Table 1. 
In the first contrast, we found no significant differences 
in growth rate between the relatively larger WR(H) off- 
spring and the smaller WR(L) offspring (F^g = 1.48, 
P = 0.23). In the second contrast, as in our earlier analy- 
ses, we did find a significant influence of density on the 
growth rate of offspring (Fijg = 11.90, P = 0.002). The 
final contrast that asked if the effect of density differed 
on WR(H) versus WR(L) fish was not significant 
(Fi,i9 = 0.01, P = 0.93). 

For age at maturity, none of the contrasts was sig- 
nificant (WR[H] vs. WR[L]: Fi_i7 = 1.03, P= 0.32 den- 
sity: Fij7 = 1.34, P = 0.26; interaction between density 
and offspring size/maternal environment: Fij7 = 0.0, 
P = 0.98). 



For the size at maturity, as was the case for age at 
maturity, none of the contrasts was significant (WR[H] 
vs. WR[L]: Fi,i7 = 0.17, P = 0.68; density: Fi,i7 = 1.16, 
P = 0.29; interaction: Fi_i7 = 1.79, P = 0.19). 



Discussion 

There were two main outcomes of this study. First, all of 
the offspring grew significantly more slowly in the higher 
density regardless of their identity and the identity of the 
competitor. Second, fish that were significantly larger at 
birth and were derived from a population that experi- 
ences extremely high densities (WR[H] fish) were better 
competitors, in the sense that, when paired with fish from 
the low-density population (TP), they caused a delay in 
the onset of sexual maturity in TP females compared with 
the onset in a homogeneous cohort of TP fish. 

Density effects in single stock treatments 

Quadrupling density decreased juvenile growth rates by 
about 50%, regardless of the population. The lack of a 
difference between populations in response to density is 
striking because fish in Trout Pond experience densities 
that are consistently two or more orders of magnitude 
lower than those at the Wacissa River (Leips and Travis 
1999); this is a chronic difference that might lead one to 
expect some density-dependent adaptation. Per capita 
resource availability in the natural populations remains 
unknown, and so it is possible that the resources available 
per individual are similar in the two populations; how- 
ever, this seems unlikely given what we know about them 
(unpubl. data from otoliths, etc.). An additional caveat is 
that some aspect of the laboratory treatment (e.g., the use 
of flake food) may have masked population-specific 
responses of growth rate to density. Future work on this 
system could test these possibilities by using more natural 
food types (e.g., zooplankton) across a wider range of per 
capita resource levels. 

There were no obvious effects of density on the size at 
sexual maturity for females. There were also no significant 
effects of density on the age at maturity for females, 
although there was evidence of a population-specific 
response to density. WR fish appeared to mature earlier 
at higher density while maturation rates of TP fish 
appeared unaffected by increased density in single stock 
treatments. We discuss these potential population-specific 
responses to density in more detail below. The lack of an 
effect of density on the age at maturity is surprising given 
its effect on early growth rates and the normally depressant 
effects of density on development rates of ectothermic 
vertebrates (Wilbur 1980). However, because we only 
measured age at maturity once a week, if the treatments 
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produced differences of less than 1 week, we would not 
have detected them. 

Equivalence of competitors? Single stock 
versus mixed stock comparisons 

With respect to initial growth rates, our results indicate 
that fish from both populations are equivalent competi- 
tors. In contrast, we found evidence of competitive asym- 
metry effects on the age at maturity, between fish from 
TP and WR. The observed patterns in age at maturity 
suggest that when placed in mixed stock treatments, fish 
from the two populations have different norms of reac- 
tion to density. If anything, WR(H) offspring appear to 
mature slightly (although not significantly) earlier when 
exposed to increased densities in every case, regardless of 
the identity of the competitor. In contrast, as density 
increased, offspring from TP delayed maturity when 
mixed with WR(H), but matured at the same age when 
in competition with members of their own population; 
we also note that there was no effect of density on growth 
rate or size at maturity of TP fish when compared across 
these treatments. So if growth rates were comparable but 
age at maturity differed, there should have been an effect 
on the size at maturity. This means that either the effects 
of WR fish on the age at maturity of TP offspring are 
overestimated by the sampling method we used, or that 
growth rates of developing females in the different treat- 
ments diverged after growth rate estimates were made 
(after the 2-week period). Another possibility is that small 
differences in growth rate may produce large differences 
in the age at maturity. There is no way to distinguish 
between these hypotheses with the current data, but these 
possibilities would be useful to explore in future experi- 
ments. 

While competition may explain the density-dependent 
depressant effects of WR fish on the age at maturity of 
TP fish, the WR fish accelerated the age at maturity with 
increasing density, regardless of the identity of the com- 
petitor. This suggests that density-dependent competition 
is not the sole factor influencing this trait in WR fish. 
One possible explanation is that density-dependent social 
effects influenced the age at maturity in WR females. 
Social effects on age at maturity are commonly reported 
in fish (Borowsky 1978; Bushman and Burns 1994; Danyl- 
chuk and Tonn 2001; Aday et al. 2003; Walling et al. 
2007) and are attributed to a variety of potential mecha- 
nisms including the existence of dominance hierarchies, 
aggressive interactions, and chemical signaling (Aday 
et al. 2003). As WR fish typically experience much higher 
densities than TP fish, it is possible that these populations 
have evolved different sensitivities to social interactions 
that influence density-dependent maturation rates. 



Adaptive Differentiation in Offspring Size 

Fitness consequences of differences in 
offspring size 

Our results indicate that under competitive conditions, 
for H. formosa, it is always better to be relatively bigger 
than your competitor. The larger offspring from the Wac- 
issa River were never worse competitors than TP offspring 
and, under some conditions, were clearly superior (TP 
delayed maturity when WR[H] were present, compared 
to TP alone and, conversely, WR[H] were unaffected by 
presence of TP). As delayed maturity increases time to 
first reproduction, all else being equal, delaying the age at 
sexual maturity should reduce fitness. Along these same 
lines, we note that the relatively larger WR(H) offspring 
matured earlier than the smaller WR(L) offspring when 
they were mixed with TP (Table 1). Although not signifi- 
cant, the trend reinforces the hypothesis that WR(H) 
offspring are better competitors than the smaller WR(L) 
offspring, and that larger offspring size may allow earlier 
maturation than smaller offspring size. In this light, we 
propose that the increased size of offspring produced by 
females experiencing high density (Leips et al. 2009, this 
study) represents an adaptive, plastic response that 
increases offspring fitness in high density, competitive 
environments. Whether the fitness advantage of being rel- 
atively larger results from competitive superiority and/or 
some other mechanism such as increased aggressiveness 
toward relatively smaller individuals remains an open 
question. 

Our results raise an important question that remains 
the subject of ongoing debate: if increased size of off- 
spring at birth is adaptive, for whom is it so (Wolf and 
Wade 2001; Marshall and UUer 2007; UUer and Pen 
2011)? More precisely, is this effect an adaptation of 
mothers or offspring? Maternal adaptation is the more 
familiar hypothesis (Marshall and UUer 2007). This argu- 
ment has two parts. First, in the context of this study, 
crowding has selected for females to allocate more 
resources to offspring because crowding signals resource 
shortage and a shortened lifespan of adults. Life history 
theory predicts that an increase in the mortality rates of 
older age-classes should select for an increase in investment 
in reproduction at ages before the increased mortality 
(Charlesworth 1994). The same logic applied to plasticity 
in female investment argues that when environmental 
cues indicate that the risk of mortality for an individual 
is about to increase, she should increase her investment 
in reproduction. Because H. formosa are matrotrophs and 
cannot quickly change their fecundity in response to 
resource level availability, we might expect the change in 
investment be reflected in offspring sizes rather than 
number (Pollux and Reznick 2011). Second, females may 
use social cues as an indicator of future high levels of 
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crowding and resource shortages, allowing them to com- 
pensate for deteriorating conditions. Isolated H. formosa 
females receiving reduced food levels have been reported 
to produce smaller offspring (Reznick et al. 1996; but see 
Travis et al. 1987). This could be maladaptive in situa- 
tions where there is selection for a large body size. There- 
fore, crowding may offer a cue about resource shortages 
that acts in advance of the actual effect of that shortage. 
This could allow females to adjust their investment 
almost immediately so that they can keep their offspring 
at a relatively constant size despite falling per capita food 
levels. This would be an example of a "counter gradient" 
phenomenon (Levins 1969; Conover and Schultz 1995); 
in nature, increased crowding would be associated with 
lower per capita food and the response we have docu- 
mented in females may act to keep offspring size near the 
optimal level. 

Increased size at birth could also/instead be an adapta- 
tion of offspring. Being larger can not only benefit 
offspring fitness, especially in competitive conditions, if 
larger size enhances resource acquisition during the juve- 
nile stage this could also enhance the reproductive fitness 
of aduks (Lindstrom 1999; Auer 2010). We know that 
offspring can influence maternal investment (Schrader 
and Travis 2009), and so we can posit that offspring can 
coerce greater maternal investment when mothers experi- 
ence higher densities. This hypothesis requires 1) that 
there be a physiological cue associated with higher 
density, 2) that the cue passes through the placenta in the 
female ovary to the embryo, and 3) that the embryo 
responds appropriately. With respect to 1), a common 
response to crowding in other vertebrates is increased 
secretion of Cortisol (RotUant and Tort 1997; Ramsay 
et al. 2009). It remains to be demonstrated whether 
female H. formosa respond similarly and whether embryos 
can seize upon this cue and increase production of com- 
pounds like IGF-2, which is associated with increased 
resource transfer from mother to embryo, to coerce 
greater investment from their mothers. If embryos do not 
respond to elevated Cortisol, females might direct 
resources away from reproduction, a commonly observed 
trend in egg laying fish (McCormick 2006; Mileva et al. 
2011). There is evidence that even exposure of unfertilized 
eggs to elevated but physiologically relevant levels of Cor- 
tisol can influence behavioral and morphological traits of 
juvenile fish (Sloman 2010) so our hypothesis is plausible. 

Our results suggest that genetically based differences in 
offspring size between populations are maintained in part 
by the selective advantage of being relatively large in a 
competitive environment (e.g., the high-density popula- 
tion of the Wacissa River). This interpretation is consis- 
tent with studies across a broad range of taxa that have 
found that larger propagule sizes are selectively favored in 



competitive environments (Hutchings 1991; Guisande 
1993; Sinervo et al. 2000; Marshall et al. 2006; Allen et al. 
2008). We recognize that additional ecological factors like 
predation may contribute to the divergence in offspring 
size between these populations. Indeed, in a 10 year study 
of life history variation in about a dozen populations of 
H. formosa (including the TP and WR populations), Sch- 
rader and Travis (2012) found that increased density and 
decreased predation risk were each independently associ- 
ated with larger offspring sizes. This result indicates that 
while density is not the only likely agent of selection on 
offspring size, density is acting independently of the 
action of predation as an agent of selection. 

A number of mechanisms can contribute to the com- 
petitive advantage of larger size; larger individuals may be 
more capable of sustained foraging activity (Parkos and 
Wahl 2010), they may be more effective at rapidly har- 
vesting scarce algal resources, be able to eat a greater 
diversity of food types, or be more likely to prevail in 
competitive interactions for resources (e.g., because of 
higher burst speeds (Segers and Taborsky 2011)). Size can 
be an important factor determining the outcome of 
aggressive encounters (e.g., Farr 1989). Whether large size 
is genetically correlated with another trait such as evolved 
differences in aggressive behavior (e.g., Herczeg and Vali- 
maki 2011) or egg composition (e.g., Ruuskanen et al. 
2011) that improves competitive ability is not known. 

While being large may also be advantageous in non- 
competitive environments, because of the tradeoff 
between offspring size and offspring number, the relative 
advantage of having large offspring in noncompetitive 
environments may not offset the resultant reduction in 
parental fitness due to the decrease in offspring number. 
This may well be the case for the TP population. Infor- 
mation on the size-specific mortality rates of juveniles in 
natural populations that differ in density could help to 
test this idea. Additional information on the relationship 
between offspring size and fitness across a range of densi- 
ties (cf Winn and Miller 1995) could also contribute to 
our understanding of the adaptive significance of patterns 
of population density and offspring size in natural popu- 
lations. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. A photo of one cattle tank with the mesh bags 
suspended in the tank. Each bag contained one replicate 
of each of the 8 treatments. 

Figure S2. Mean growth rates of offspring (± 1 SE) from 
single stock treatments (means for treatment groups com- 
pared in Analysis A). TP = Offspring of Trout Pond 
Females held at high maternal density, WR(H) = Offspring 
of Wacissa River females held at high maternal density. 
Figure S3. Mean age at maturity of offspring (± 1 SE) 
from single stock treatments (means for treatment groups 
compared in Analysis A). TP = Offspring of Trout Pond 
Females held at high maternal density, WR(H) = Offspring 
of Wacissa River females held at high maternal density. 
Figure S4. Mean size at maturity of offspring (± 1 SE) 
from single stock treatments (means for treatment groups 
compared in Analysis A). TP = Offspring of Trout Pond 
Females held at high maternal density, WR(H) = Offspring 
of Wacissa River females held at high maternal density. 
Figure S5. Mean growth rates of TP offspring (± I SE) 
single stock versus mixed stock treatments (means for 
treatment groups compared in Analysis B). TP = Offspring 
of Trout Pond Females held at high maternal density, WR 
(H) = Offspring of Wacissa River females held at high 
maternal density. 

Figure S6. Mean growth rates of WR(H) offspring (± 1 
SE) single versus mixed stock treatments (means for treat- 
ment groups compared in Analysis C). TP = Offspring of 
Trout Pond Females held at high maternal density, WR 
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(H) = Offspring of Wacissa River females held at high 
maternal density. 

Figure S7. Mean size at maturity of WR(H) offspring 
(± 1 SE) single versus mixed stock treatments (means for 
treatment groups compared in Analysis C). TP = Offspring 
of Trout Pond Females held at high maternal density, WR 
(H) = Offspring of Wacissa River females held at high 
maternal density. 

Table SI. The top four rows reexpress the information in 
Table Al. The top three rows indicate the identity and 
density of fish in that treatment level (Population 
MD = population and maternal density) and a fourth 
row that indicates the cell number used in Table Al for 
that combination of TP (H), WR(H), or WR (L) fish. 
The remainder of the table is a set of contrasts among 
treatment levels that describe four two-way analyses of 
variance with different combinations of levels used in 
each analysis. The four analyses are designated by the 
Roman letters A-D. Within any single "analysis" (A, B, 
C, or D), there are three contrasts designated i, ii, and iii. 
The first two contrasts, i and ii, test for main effects and 
the third contrast, iii, tests for the interaction between 
two main effects. Analyses and hypotheses tested contrasts 
within each analysis are described in the text. 
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Table Al. The experimental matrix of treatment combinations. Num- 
bers in each cell refer to designations for contrasts between different 
treatments. Densities: 1/2A/ = 4 individuals, W = 8 individuals and 
4A/ = 32 individuals. T = Trout Pond offspring, W = Wacissa River off- 
spring (H = offspring of females maintained at high density, L = off- 
spring from females maintained at low density). 
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